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SYSTEM  TO  MEASURE  CARRIER  DRIFT  VELOCITY  IN  SEMICONDUCTOR  MATERIALS 


1.  INTRODUCTION 


The  electrical  behavior  of  semiconductor  devices  is  strongly 
dominated  by  the  motion  of  charge  carriers  in  the  electric  field 
present  within  the  device.  Particle  motion  in  a  crystalline  solid 
differs  from  that  in  free  space,  for  example,  because  the  energy 
transferred  from  the  field  to  the  particle  is  influenced  by  the 
energy  band  structure  of  the  material,  scattering  processes  which 
transfer  energy  from  the  particle  to  the  crystal  lattice  and  a 
variety  of  effects  such  as  trapping,  intra-  and  inter- valley 
scattering,  etc.  One  means  of  accounting  for  these  complex  effects 
is  to  define  an  average  drift  velocity  for  carriers  in  the  presence 
of  an  electric  field.  For  some  III-V  semiconductors,  such  as  GaAs, 
the  velocity-field  curve  is  non-monotonic  and  exhibits  a  region  of 
negative  differential  mobility  as  shown  in  Figure  1.1. 

The  negative  differential  mobility  region  on  the  high  field 
side  of  the  peak  gives  rise  to  unique  phenomena  such  as  Gunn 
oscillations  (1)*.  Precise  knowledge  of  the  drift  velocity-electric 

*  The  numbers  in  the  parentheses  in  the  text  indicate  reference  in 
the  Bibliography. 
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1  Reported  velocity-field  characteristic  for  GaAs  (3) 
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field  relationship  is  needed  to  understand  device  behavior  and  tu 
aid  in  the  design  and  improvement  of  semiconductor  devices  such  as 
field  effect  transistors  and  Gunn  diodes. 

The  problem  of  measuring  the  velocity-field  curve  can  be 
illustrated  by  considering  the  behavior  of  charge  in  a  vacuum  diode. 
Figure  1.2  shows  a  parallel  plate  diode  whose  plates  are  separated 
by  a  length  L  and  with  a  constant  internal  electric  field  equal  to 
the  applied  voltage  divided  by  the  length.  If  a  charge  layer  of 
electrons,  a  is  injected  between  the  plates  it  will  induce  a 
surface  charge  c^at  x=0  and  o2  on  the  contact  at  x=L  where  oi  and 
<J2  are  given  by 


L~*l 

aI  *  -°s  T 


■  -<?  -r4- 

2  s  L 


(1.1) 

(1.2) 


With  the  applied  voltage  shown  the  negative  charge  sheet,  a  , 
will  move  to  the  right  with  some  velocity  v  and  the  charge  densities 
at  the  terminals  will  change  with  time,  inducing  current  in  the 
external  circuit  related  to  the  current  density  JT: 


d_!i  jl., 

dt  dt  (-as  L  ^ 


(1.3) 
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Here  p  is  the  volume  charge  density  and  v  is  the  free  space 
velocity.  For  most  experiments  in  determining  v(E)  this  external 
circuit  current  is  the  only  observable  variable. 

In  a  vacuum  diode  with  a  uniform  electric  field  the  velocity  v 
increases  with  time.  However,  in  a  semiconductor  the  velocity  will 
depend  on  the  v( E)  relation.  If  a  uniform  field  is  present  the 
velocity  will  be  constant.  In  a  typical  experiment  charge  injected 
into  the  diode  region  induces  a  pulse  of  current  in  the  external 
circuit.  The  pulse  rise  time  is  proportional  to  the  rate  of  charge 
injection.  The  pulse  width  measures  the  transit  time  of  charge 
through  the  diode.  The  pulse  fall  time  is  related  to  trapping  and 
diffusion  of  the  charge  layer  during  transit. 

To  measure  the  velocity  within  a  semiconductor  medium,  a 
nonconducting  Schottky  contact  is  usually  employed.  A  metal,  such 
as  A1 ,  is  evaporated  on  an  epitaxial  semiconductor  layer  grown  on  a 
conducting  substrate  that  forms  the  anode.  The  difference  between 
the  metal  work  function  and  the  n-type  semiconductor  work  function 
will  cause  electrons  to  travel  from  the  material  into  the  metal. 
This  leaves  a  "depleted"  region  beneath  the  metal  which  has  a  net 
positive  charge  and  a  negative  charge  on  the  metal.  An  electric 
field  is  created  with  its  maximum  (negative)  value  at  the  metal.  It 
increases  to  zero  within  the  semiconductor.  The  electron  migration 
stops  when  the  electrostatic  potential  of  the  semiconductor  is 
raised  to  match  the  work  function  differential  across  the  interface. 
Under  increasing  reverse  bias  this  depletion  region  grows  into  the 
material . 


For  a  reverse  biased  Schottky  diode  the  situation  is  similar  to 
the  vacuum  diode  case  except  that  the  cathode  is  now  the  Schottky 
metallization,  the  anode  is  the  end  of  the  depletion  region  and  the 
electric  field  slope  is  related  to  the  local  doping  density,  Nd(x), 
by  Poisson's  equation 


dE  =  £ 
dx  "  £ 


Nd(x) 


(1.5) 


where  q  is  the  electron  charge  and  c  the  dielectric  permittivity  of 
the  semiconductor  medium.  The  electric  field  will  therefore  go  from 
a  maximum  intensity,  Eq,  at  the  edge  of  the  Schottky  contact  to  zero 
past  the  end  of  the  depletion  region.  Figure  1.3a  shows  the  field 
profile  for  an  arbitrarily  doped  Schottky  diode  of  length  L  under 
normal  bias.  Some  epitaxial  layers  are  short  enough  and 
sufficiently  lightly  doped  that  the  electric  field  extends  into  the 
substrate  which  has  a  much  higher  doping  concentration.  In  the 
substrate  the  field  drops  off  very  quickly  with  depth  making  a 
profile  which  looks  similar  to  Figure  1.3b.  This  is  called  the 
punched- through  case  where  the  electric  field  has  a  non-zero  value 
at  the  end  of  the  epitaxial  layer.  The  contribution  of  the  carrier 
transit  time  in  the  substrate,  in  this  case,  is  negligible  compared 
to  that  of  the  epitaxial  layer  since  the  length  of  field  penetration 
is  minute  in  the  heavier  doped  material.  Most  velocity-field 
measurements  are  done  with  the  sample  biased  in  the  punched- through 


state. 


f(x) 


vf(£)  H 


3  Electric  field  profile  (a) 
(b)  punched  through 


Earlier  techniques  to  determine  velocity-field  relationships 
used  Schottky  diodes  with  long,  lightly  doped  material.  A  pulse  of 
electrons  was  injected  at  one  end  with  an  electron  gun  and  the 
transit  time  of  the  charge  across  the  depletion  region  was  measured. 
This  was  known  as  the  time-of-fl ight  technique  (2). 

Later  a  modulated  sinusoidal  signal  injection  scheme  was 
developed  to  yield  better  resolution  (3, 4, 5, 6).  The  phase  of  the 
detected  signal  contains  information  concerning  the  transport  of 
charge  particles  through  the  lattice.  If  sufficiently  high 
modulation  frequencies  are  used  the  semiconductor  sample  could  be 
small  compared  to  those  used  for  time-of-fl ight  measurements  and 
thus  more  nearly  approximate  conditions  found  in  actual  devices. 
These  microwave  phase  (CW)  time-of-fl ight  schemes  have  been  used 
before  and  show  much  promise  for  small  samples  with  high  doping 
levels. 

Measuring  the  time-of-fl ight  mobility  of  semiconductors  is 
difficult  and  requires  a  large  amount  of  equipment.  A  means  of 
injecting  free  carriers  at  a  constant  rate  and  known  modulation  into 
the  active  region  is  needed.  Target  samples  of  proper  doping 
profile  and  epitaxial  layer  thickness  are  necessary  for  detailed 
analysis  of  sections  of  the  v(E)  curve.  Sensitive  signal  detection 
circuitry  is  required  along  with  an  accurate  technique  of 
interpreting  the  measured  phase-voltage  data  and  deriving  the 
velocity-field  relationship. 
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Attempts  have  been  made  in  the  past  to  determine  the  mobility 
of  certain  semiconductors.  Some  of  the  previously  stated 
requirements  have  been  considered. 

1.1  PREVIOUS  MEASUREMENT  APPROACHES 


Electron  beam  injection  methods  were  employed  by  J.G.  Ruch  and 
G.S.  Kino  (2)  in  the  time-of-fl ight  measurement  scheme.  In  their 
experiment  the  beam  was  swept  across  an  aperture  exposing  the  target 
material  and  the  duration  of  the  current  pulse  in  the  external 
circuit  divided  by  the  sample  length  was  taken  as  the  velocity. 
Sample  doping  was  compensated  low  so  a  uniform  electric  field 
assumption  could  be  made  and  long  samples  (300  microns)  were  used 
to  get  measurable  pulse  widths.  These  long  samples  allowed 
diffusion  effects  to  become  significant,  enabling  diffusion  vs. 
electric  field  measurements  but  also  increasing  the  pulse  fall  time. 
Uncertainty  in  the  pulse  duration  along  with  the  lack  of  system 
calibration  encouraged  the  investigation  of  other  measurement 
techniques. 

In  1974  Evans  and  Robson  (3)  reported  data  taken  using  a 
microwave  phase  time-of-fl ight  method.  Figure  1.4  shows  the 
schematic  diagram  of  their  experimental  setup.  The  beam  was 
deflected  past  the  sample  face  at  an  RF  rate  and  the  phase  shift  of 
the  external  circuit  current  was  measured.  The  doping 
concentrations  were  low  and  the  sample  lengths  short  so  a  uniform 


Figure  1.4  Evans  and  Robson's  microwave 
time-of-f light  apparatus  (3). 
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electric  field  equal  to  the  voltage  divided  by  the  epitaxial  region 
length  was  assumed.  Since  the  electric  field  was  taken  to  be 
constant  so  was  the  velocity  which  was  calculated  from  the  measured 
phase  using  the  expression 


where  w  is  the  microwave  radian  frequency  of  the  current.  This 
result  will  be  derived  later. 

The  CW  method  and  uniform  field  approximation  are  being  used  at 
both  Cornell  University  (6)  and  the  University  of  Illinois  (4,5)  on 
Si,  GaAs  and  InGaAs.  Besides  using  low  doped  samples  most  of  these 
measurements  are  being  performed  at  high  electric  fields  where  the 
velocity  is  more  nearly  constant  to  minimize  velocity  averaging 
errors. 

Monte  Carlo  simulations  (7)  using  band  structure  probabilities 
have  predicted  the  doping  and  temperature  dependence  of  the  v(E) 
curve  for  GaAs  and  InP.  These  results  can  be  compared  to  measured 
data  at  high  doping  levels. 

The  work  here  deals  with  the  measurement  of  the  v(E)  curve  at 
low  electric  fields.  It  addresses  and  analyzes  the  non-uniform 
field  situation,  signal  detection  considerations  and  numerical 
methods  to  relate  the  measured  current-phase  to  the  velocity- field 
relationship. 
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1.2  SCOPE  OF  THIS  WORK 


The  system  employed  in  the  work  reported  here  is  an  improved 
model  of  the  one  described  in  ONR  report  #80-2  (8)  which  is  similar 
in  design  to  that  of  Evans  and  Robson  (3).  Figure  1.5  shows  a  block 
diagram  of  the  measurement  system.  An  electron  beam  is  modulated  at 
6  GHz  and  swept  past  a  Schottky  diode  face.  Carriers  are  injected 
through  the  Schottky  contact  into  the  sample  with  a  frequency  of  12 
GHz  since  the  beam  passes  the  sample  aperture  twice  in  one 
modulation  period.  The  phase  of  the  resulting  signal  in  the 
external  circuit  is  then  measured  as  a  function  of  applied  reverse 
bias,  relative  to  the  reference  signal  phase,  with  a  Hewlett-Packard 
automatic  network  analyzer.  Chapter  3  contains  a  detailed 
description  of  the  system  components  and  measurement  technique. 
There  are  a  number  of  advantages  this  analysis  has  over  the 
previously  mentioned  microwave  time-of-flight  measurements. 

The  doping  profile  of  the  material  used  in  these  measurements 
is  well  known  and  so  the  electric  field  can  be  calculated  with 
respect  to  depth  and  applied  bias.  A  nonuniform  field  solution  is 
possible  using  numerical  methods  to  derive  a  v(E)  curve  which  best 
fits  the  observed  data.  This  technique  also  allows  the  examination 
of  higher  doping  concentrations  than  could  be  done  using  previous 
techniques. 

This  scheme  should  yield  information  about  the  v(E)  curve  down 
to  zero  electric  field.  These  measurements  are  then  absolute  and  no 


1.5  Block  diagram  of  v(E) 
experimental  apparatus. 
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other  data  are  needed  as  reference.  A  Hewlett-Packard  ANA  under 
computer  control  is  used  to  automate  the  measurements. 

The  remainder  of  this  report  is  devoted  to  describing  the 
theorectical  aspects  of  the  measurement,  the  physical  design  of  the 
system,  and  the  data  generated  using  this  technique.  Results  are 
compared  to  previous  experiments  and  simulations.  Conclusions  and 
proposals  for  improvements  are  made  from  these  comparisons. 
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2.  THEORETICAL  BASIS  FOR  V(E)  MEASUREMENTS 


Through  measurement  and  interpretation  of  the  external  circuit 
current  the  carrier  velocity- internal  electric  field  relation  can  be 
determined.  The  phase  of  this  current  holds  enough  information  to 
uniquely  define  the  v(E)  curve. 

2.1  THE  CARRIER  VELOCITY  -  CURRENT  PHASE  RELATIONSHIP 


First  the  current  continuity  equation  within  the  material  must 
be  examined 


(2.1) 


This  can  be  simplified  if  only  current  density  variations  in  the 
direction  of  the  applied  electric  field,  say  along  x,  are  considered 


M  _  i£_ 

ax  at 


(2.2) 


The  range  of  x  is  from  the  Schottky  metal ization  (x=0)  to  the  end  of 
the  depletion  region. 


The  electron  beam  is  swept  past  the  sample  aperture 
periodically  at  a  frequency  of  w/2  as  described  in  Chapter  1.  The 
free  carrier  density  and  the  current  density  can  be  assumed  to  vary 
as  exp(jwt)  since  there  are  two  current  pulses  per  modulation  cycle. 
Equation  2.2  can  thus  be  written  as 

||+j-o-0  (2.3) 


Substituting 


'V. 

J 


'V 

*  pv 


(1.4) 


into  (2.3)  yields 


a3r  } 

—  +  iu>— 
3X  J  v 


(2.4) 


a  first  order  differential  equation  relating  the  local  current 
density  and  the  local  drift  velocity,  neglecting  diffusion  effects. 
The  homogeneous  solution  is 

J-J0«xP  <2-5> 


with  J0  as  a  constant  dependent  on  the  number  of  injected  electrons 
and  their  energy.  Using  Poisson's  equation 
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and  substituting  into  the  current  continuity  equation 


(2.7) 


the  total  current  can  be  defined  as 


J 


terminal 


J  +  e 


jiE 
d  t 


(2.8) 


which,  if  integrated  over  the  active  region  of  the  diode,  relates 
the  drift  velocity  to  the  measurable  external  circuit  current.  The 
applied  voltage  is  assumed  to  be  constant  so  the  electric  field  can 
be  taken  to  be  time  invariant.  Thus  the  terminal  current  is  given 
by 


(2.9) 


The  preceeding  derivation  was  performed  by  Evans  and  Robson  (3) 

to  relate  terminal  current  phase  to  localized  drift  velocity.  To 

further  simplify  the  measurement,  samples  with  low  doping 

concentrations  can  be  used  so  the  electric  field  and  the  velocity 

can  be  considered  uniform  along  the  sample.  With  a  constant 
velocity  (2.9)  becomes 

i 

exp  (  -^^)dx 

V 


I 


(2.10) 


which  can  be  expressed  as 


% 
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h.  i.  (exp  -hJ. 
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(2.12) 


The  phase  of  the  external  circuit  current  is 


<e 


2v 


(1.6) 


and  the  amplitude  of  the  alternating  current  density  is 


A 


!  J  I  Isinc  ~ 
o  2v 


(2.13) 


When  a  bias  voltage  is  applied  to  a  sample  under  the  uniform  field 
assumption  the  internal  electric  field  is  equal  to  the  voltage 
divided  by  the  depletion  length.  From  the  measured  current  phase 
the  velocity  can  be  calculated  using  (1.6).  Thus  points  are 
generated  on  the  v(E)  curve. 
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2.2  NONUNIFORM  ELECTRIC  FIELD  SOLUTION 


For  samples  with  doping  concentrations  >5x10  the  field  slope 
will  be  greater  than  6.8  kV/cm  per  micron  so  for  a  3  or  4  micron 
depletion  region  the  change  in  electric  field  can  be  greater  than  20 
kV/cm.  The  uniform  field  approximation  can  not  be  used  except  for 
areas  of  low  mobility,  where  the  velocity  is  near  constant  with 
respect  to  electric  field. 

For  a  given  doping  profile  and  applied  voltage  the  electric 
field  profile  of  the  sample  can  be  calculated.  For  a  monotonical ly 
decreasing  electric  field  the  external  circuit  current  phase  can  be 
expressed  as  an  integral  over  an  electric  field  interval  on  the.  v(E) 
curve  instead  of  integrating  over  the  depletion  length.  An  initial 
guess  about  the  v(E)  curve  can  be  made  and  corrected  using  iterative 
numerical  methods  to  yield  a  v(E)  relationship  which  satisfies  the 
measured  current  phase.  To  derive  such  a  phase-v(E)  relationship 
(2.9)  must  be  restated  in  terms  of  the  terminal  current  phase  as 


<P  -  Arg 


i  r  r * 

)  f  [  f  dx  1 

j  J  exppuJ  ^00  j 


(2.14) 


where  is  the  epitaxial  layer  depth  in  the  punched- through  case  or 
the  point  where  the  electric  field  goes  to  zero  in  non-punched 
through  situations.  See  Figure  1.3.  From  examining  the  doping 
profile  (2.14)  can  be  changed  from  an  integration  over  length  to  an 
integration  over  electric  field. 


L*. 


I 
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relationship.  The  key  to  understanding  the  analysis  method  used 

here  is  to  recognize  that  each  measurement  describes  a  double 

integration  on  a  specific  interval  of  the  v(E)  curve. 

As  the  electron  passes  through  the  depletion  region  it  moves 

with  a  velocity  related  to  the  local  electric  field.  For  a  given 

applied  voltage  the  electric  field  can  be  determined,  along  with  its 

most  negative  value  at  the  metallization,  and  its  least  negative  at 

the  end  of  the  epitaxial  layer,  assuming  a  punched- through  state. 

Figure  2.1  illustrates  a  typical  v(E)  curve  divided  into  measurement 

intervals,  each  with  a  corresponding  current  phase.  The  external 

current  phase  is  a  double  integration  on  each  section  of  the  curve 

so  a  method  must  be  used  to  derive  the  curve  which  best  satisfies 

the  measured  phases.  To  perform  this  de-embedding,  electric  field 

intervals  can  be  chosen  to  yield  sufficient  resolution  around  the 

_  v(E)  curve  regions  of  greatest  mobility  deviation.  The  integration 

i 

of  the  curve  masks  information  about  its  shape  within  the  electric 
field  limits.  However,  many  integrations  over  different  intervals 
can  give  enough  information  so  that  a  correct  curve  can  be  obtained. 

From  lumped  element  modeling,  which  will  be  explained  in 
Chapter  3,  punched-through  biasing  was  found  to  be  the  only  way  to 
perform  these  measurements.  Samples  with  epitaxial  layer  doping 

15  -3 

concentrations  on  the  order  of  2x10  cm  and  3  to  5  microns  in 
length  were  necessary  for  good  analysis  of  the  raw  data.  The 
program  used  for  the  analysis  is  listed  in  Appendix  7.1  and  is 
described  in  block  diagram  form  in  Figure  2.2.  A  general 
description  of  its  function  is  as  follows: 


ME  AS  * 2 

Electric  Field 

Figure  2.1  Electric  field  intervals  for  individual 
phase  measurements 
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Figure  2.2  Block  diagram  of  Newton-Raphson 
iterative  method. 
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Figure  2.2  (cont.) 


output  v(E) 


Figure  2.2  (cont.) 


a)  An  initial  guess  is  made  about  the  shape  of  v(E)  curve  and 


numerical  integration  is  performed  to  produce  the  expression  G(e) 


(2.19) 


b)  Euler’s  relation  allows  equation  (2.18)  to  be  represented  by 
two  integratable  expressions,  one  real  and  one  imaginary. 

$  *  Argj  I  cos(uG(c))h'(e)dc  -j  sin(wG(*))  h’(e)d£  (220) 

E0 

c)  After  integration  the  half  angle  relation  is  used  to  find  the 
calculated  phase  and  a  percentage  error,  L,  is  found  between  the 
measured  and  the  calculated  phase  values. 


L  a  $meas_2_£cale 
$meas 


(2.21) 


d)  A  measurement  electric  field  interval  on  the  v(E)  curve  is 
incremented  by  a  constant  velocity  step  across  the  range.  Steps 
a)  through  c)  are  then  repeated  and  the  change  in  percent  phase 
error  for  the  velocity  increment  can  be  found  for  each  measured 
point.  This  interval  incrementing  is  done  for  each  measurement 
with  the  percent  error  changes  calculated  for  all  the  phases  each 
time  until  a  Jacobian  matrix,  J,  is  generated. 


B  . -g^ . —  — -  . 
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The  velocity-field  curve  is  then  augmented  by 

-1 

v*  ■  v  -  f  *  J  (2.22) 

where  the  I  matrix  is  the  percentage  difference  between  the 
measured  phase  and  the  phase  calculated  from  the  previous  v(E) 
curve. 

e)  Steps  a)  through  d)  are  repeated,  iterating  the  v(E)  guess  each 
time,  until  the  percentage  error  for  each  measured  point  falls 
below  some  minimum  value.  The  resulting  velocity-field  curve  is 
then  the  unique  solution  for  the  measured  data. 

The  above  method  assumes  zero  diffusion,  no  phase  ambiguity, 
and  accurate  phase  measurements.  Diffusion  was  neglected  primarily 
because  the  analysis  becomes  too  difficult  otherwise.  The  sample 
lengths  were  kept  short  to  minimize  diffusion  effects  on  the  phase 
of  the  detected  signal.  This  method  can  be  used  for  any  v(E) 
relation  and  any  semiconductor  where  diffusion  is  not  too  great 
(<15%  phase  innacuracy). 

The  accuracy  of  the  raw  phase-voltage  data  limits  the  accuracy 
of  any  measurement  technique.  Therefore  a  system  had  to  be  designed 
to  generate  reliable  raw  data.  Chapter  3  describes  the  hardware 
used  in  this  scheme  to  produce  useful  phase  vs.  voltage 


measurements. 
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3. SYSTEM  REALIZATION 


Performing  the  phase  vs.  applied  voltage  measurements  requires 
a  dependable  carrier  injection  and  current  phase  detection  system. 
Each  part  of  this  system  had  to  be  fabricated  and  individually 
tested  to  meet  our  specific  criteria  of  performance  and 
reproducability.  A  detailed  description  of  the  electron  gun,  power 
supply  and  second  stage  amplifier  is  available  to  the  reader  in 
O.N.R.  Report  #80-2  (8).  Machine  drawings  for  certain  sub-systems 
and  circuit  diagrams  for  the  power  supply  can  be  found  here  in  the 
appendices.  Figure  1.5  is  a  block  diagram  of  the  entire  apparatus. 
The  following  sections  describe  the  operation  of  the  individual 
sub-systems. 

3.1  ELECTRON  GUN  ANO  POWER  SUPPLY 


A  means  of  injecting  a  large  number  of  free  carriers  into  the 
depletion  region  of  the  Schottky  diode  had  to  be  designed  with  high 
stability  and  minimal  unwanted  modulation.  Electron  beam  injection 
was  chosen  since  the  technology  is  well  known  and  the  technique  is 
applicable  to  almost  all  semiconductor  materials.  Constant  beam 
position  with  varying  beam  potential  or  beam  penetration  into  the 
sample  is  also  desired. 
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Figure  3.1  shows  a  sketch  of  the  HP  500  MHz  oscilliscope  gun 
used  as  the  electron  source.  The  potentials  relative  to  the  mesh 
cup  are  indicated.  The  filament  is  a  thoriated  tungsten  wire 
mounted  on  a  ceramic  standoff  within  the  grid  cup.  The  indicated 
potentials  were  chosen  to  give  the  highest  beam  current  density  at 
the  sample  target. 

In  anticipation  of  making  measurements  at  two  different 
frequencies  the  vertical  deflection  plates  were  used  as  a  second 
modulator.  These  plates  are  actually  a  helical  ribbon  pair  designed 
to  have  a  300  ohm  characteristic  impedance  up  to  500  MHz.  In  order 
to  drive  the  vertical  modulator  a  300  ohm  transmission  line  was 
implemented  using  two  nickel  wires  separated  by  teflon  spacers. 
Modulation  power  was  coupled  to  the  vertical  pins  on  the  gun 
connector  through  a  bifiliar  wrapped  toroid  which  provides  isolation 
and  impedance  matching.  Modulation  at  500  MHz  was  observed  on  a 
phosphor  target  as  shown  in  Figure  3.2  and  a  1  GHz  signal  was 
detected  through  a  sample  diode. 

The  filament  of  the  gun  is  held  at  a  voltage  much  lower  than 
that  of  the  sample  target  which  is  at  earth  ground.  The  special 
negative  high  voltage  power  supply  shown  in  Figure  3.3  had  to  be 
built  for  this  experiment. 

Two  supplies  were  used  to  give  a  constant  voltage  difference 
across  the  zener  diode  and  resistor  chain  used  to  obtain  all  the 
required  gun  potentials.  This  allowed  the  pick  off  points  to  be 
maintained  at  the  same  voltage  relative  to  the  reference  for  a 
variable  beam  potential.  High  voltage  supply  #1  is  an  adjustable 


BE  A  M  PATTERN 
NO  MODULATION 


Figure  3.3  Hock  diagraa  of  high  voltage  auppllea 


reference  using  negative  feedback  for  regulation.  The  output  of  the 
flyback  transformer  labeled  HV#1  is  fedback  through  a  500  Megohm 
series  resistance  to  the  negative  input  of  a  741  operational 
amplifier  which  regulates  the  amplitude  of  the  current  spikes  from 
the  switching  circuit  through  the  transformer.  The  500  Megohm 
resistance  and  17.6  nF  capacitance  in  shunt  with  the  flybacks 
provide  sufficient  smoothing  of  the  high  voltage  waveform. 

The  second  high  voltage  supply  is  identical  to  the  first  except 
its  input  maintains  HV#2  at  a  fixed  potential  offset  from  HV#1, 
usually  by  -2.8  kV.  Both  transformer  switching  circuits  run  off  the 
same  asynchronous  15  kHz  signal  to  eliminate  uncoordinated  ripple  in 
the  high  voltage  output.  The  flybacks  and  smoothing  circuits  are 
housed  in  aluminum  boxes  and  coaxial  switching  circuit  lines  are 
used  to  minimize  electromagnetic  pickup  on  the  gun  controls.  The 
entire  supply  is  floating  so  a  high  value  resistor  placed  between 
floating  ground  and  earth  ground  can  be  used  to  indicate  the  beam 
current  passing  through  the  final  beam  aperture. 

The  filament  supply  circuit  board  is  similar  in  layout  to  the 
switching  circuit  board  used  for  the  high  voltage  supplies  and  is 
also  timed  with  the  HV  switching  waveform.  The  spiky  current  passes 
through  the  primary  winding  of  an  isolation  transformer,  which  is 
also  electrically  shielded,  whose  secondary  is  within  100  V  of  HV#2 
potential.  The  secondary  signal  is  then  fullwave  rectified  and 
smoothed  to  provide  two  amps  DC  current  to  the  filament  which  has  a 
cold  resistance  of  about  one  ohm.  There  is  virtually  no  observable 
density  modulation  on  the  electron  beam. 
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The  pickoff  chain  as  seen  in  Figure  3.3  was  chosen  for  ease  of 
operation  and  optimum  DC  beam  current  through  the  sample.  With  one 
plate  of  both  horizontal  and  vertical  deflection  tied  to  reference  a 
plus  or  minus  seven  volt  swing  on  the  opposite  plate  gave  adequate 
deflection  with  the  right  amount  of  sensitivity  on  the  control  knob. 
Astigmatism  and  focus  were  chosen  through  observation  for  best  beam 
geometry  and  the  filament  and  grid  potentials  were  picked  so  the 
grid  had  enough  negative  offset  to  shut  off  the  electron  beam 
completely. 

Figure  1.5  shows  the  electron  gun  and  power  supply 
configuration  with  the  rest  of  the  system.  The  gun  is  placed  in  a 
pyrex  tube  which  is  painted  both  inside  and  out  with  conducting 
paint  (an  Aquadag  E  and  sodium  silicate  mixture)  to  prevent  static 
charge  bui'»d-up  on  the  tube  surfaces.  The  interior  coating  is  tied 
to  HV  reference  while  the  exterior  is  at  earth  ground.  The  gun  and 
all  related  surfaces  exposed  to  hard  vacuum  were  cleaned  with  TCE, 
acetone  and  isopropyl  alcohol  before  pumping.  The  pyrex  tube  forms 
the  vacuum  jacket  of  the  apparatus. 

The  resulting  electron  beam  is  350  nanoamps  with  a  diameter  on 
the  order  of  30  mils.  The  range  of  beam  potential  is  3kV  to  18kV, 
but  the  usual  operating  point  is  6.75kV.  The  500  MHz  modulator 
provides  half  an  inch  deflection  on  a  phosphor  target  with  half  a 
watt  incident  on  the  bifiliar  wrapped  transformer  of  the  vertical 
deflection  input. 


3.2  6  GHz  MODULATOR 


High  frequency  modulation  of  the  charge  carrier  density  in  the 
active  region  of  the  sample  is  required  to  get  enough  phase 
resolution  to  derive  the  velocity-field  relationship.  Calculations 
have  shown  that  for  a  diode  with  a  4  urn  epitaxial  layer,  12  GHz  will 
provide  sufficient  phase  difference  for  the  HP  ANA  to  produce 
accurate  measurements  in  the  electric  field  regions  of  interest  (8). 
A  6  GHz  modulator  was  therefore  required  which  would  produce 
adequate  transverse  deflection  at  the  sample. 

For  convenience  and  effective  operation  a  reentrant  coaxial 
cavity  modulator  as  shown  in  Figure  3.4  was  designed.  The  modulator 
body  element  consists  of  a  cylindrical  bore  in  its  center  to  form 
the  outer  wall  of  the  resonator.  A  conducting  center  post  is 
located  coaxially  within  the  bore.  This  configuration  looks  like  a 
section  of  coaxial  transmission  line  near  a  quarter  of  a  wavelength 
at  6  GHz.  Figure  3.5  illustrates  the  radial  standing  wave  pattern 
within  the  coax.  The  ported  end  at  the  bottom  appears  as  a  low 
impedance  and  the  beam  aperture  appears  to  be  capacitively  loaded, 
close  to  an  open  circuit.  This  cavity  resonates  at  6  GHz  and  the 
electric  field  between  the  post  and  the  top  of  the  cavity  is  nearly 
vertical  as  shown  in  Figure  3.5.  Coupling  to  the  resonator  is 
provided  b'  an  iris  at  the  end  of  a  section  of  ridged  waveguide  as 
seen  in  Figure  3.4. 


Figure  3.4  Reentrant  coaxial  cavity  modulator  body 
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The  modulator  snaps  into  place  in  the  mesh  cup  and  is  held  at 
HV  reference.  Its  edges  are  rounded  to  minimize  arcing  and  the 
outer  conductor  walls  of  the  coax  at  the  aperture  were  narrowed  to 
increase  the  angle  of  acceptance  of  the  electron  beam.  Notches  were 
made  in  the  modulator  perimeter  to  improve  vacuum  pumping  and  the 
body  was  gold  plated  to  reduce  secondary  emission. 

Coupling  to  the  modulator  was  done  through  the  pyrex  tube  wall 
with  a  ridged  waveguide  horn.  Figure  3.6  shows  the  ridged  waveguide 
launcher  used.  A  50  ohm  coaxial  probe  attached  to  the  ridge  0.25 
inches  long  (a  quarter  wavelength  at  6GHz)  coupled  the  6  GHz  signal 
to  the  waveguide  from  coaxial  input.  The  charateristic  impedance  of 
the  guide  is  about  140  ohms  so  a  brass  slab  air  stripline  center 
conductor  was  used  to  make  a  90  ohm,  quarter  wavelength, 
transformer  to  improve  the  match  between  the  50  ohm  coax  and  the 
ridged  waveguide. 

The  mouth  of  the  launcher  and  the  modulator  were  flared  to 
improve  power  transmission  through  the  pyrex.  Both  pieces  were 
machined  to  mate  flush  with  the  tube  wall. 

At  resonance,  the  return  loss  was  measured  at  -16  dB  at  the 
launcher  input.  The  6  GHz  deflection  system  provides  from  a  tenth 
to  a  fifth  of  an  inch  deflection  at  the  sample  with  half  a  watt 
incident  on  the  launcher,  as  can  be  seen  in  Figure  3.2.  Power  is 
provided  by  a  6  GHz  Gunn  oscillator  with  an  attenuator  to  maximize 
the  12  GHz  component  in  the  detected  sample  signal  and  an  isolator 
to  minimize  frequency  drift  from  load  variation.  Detectable  12  GHz 
signals  were  obtained  with  modulation  power  as  low  as  5  mW. 
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3.3  SAMPLE  AND  SAMPLE  HOLDER  ASSEMBLY 


The  target  sample  used  in  this  system  is  a  reverse  biased 
Schottky  diode.  Figure  3.7  shows  a  sample  cross  section.  On  a 
highly  doped  substrate  approximately  100  microns  thick  a  2  micron 
buffer  with  doping  concentration  on  the  order  10  cm  is  grown 
to  prevent  back  diffusion  of  dopant  during  epitaxial  growth.  An 
epitaxial  layer  is  then  grown  between  2  and  5  microns  thick  with 
doping  in  the  range  of  4xl014  and  2xl015  cm'3. 

The  Schottky  barrier  consists  of  a  0.025  inch  diameter  ,  200  8 
thick  circle  of  aluminum  evaporated  on  the  epitaxial  surface.  To 
penetrate  this  metallization  requires  approximately  one  keV  of  beam 
energy  (7)  so  the  energy  of  the  beam  within  the  material  still  far 
exceeds  the  band  gap  energy  of  the  semiconductor.  This  excess 
energy  is  expended  in  the  generation  of  electron-hole  pairs  within 
the  first  half  micron  of  the  depletion  region.  Current 
multiplication  on  the  order  of  a  thousand  has  been  observed, 
boosting  the  detected  12  GHz  signal  by  approximately  30  dB.  A  0.005 
inch  wide,  7000  ft  thick  guard  ring  was  also  evaporated  on  top  of  the 
circle  for  bond  wire  contact  to  the  conventional  C3  package  used  to 
mount  the  sample  in  the  measurement  system.  On  later  samples  a  bond 
pad,  0.002  inches  by  0.002  inches,  was  used  in  place  of  the  guard 
ring  as  the  bonding  area.  More  current  was  detected  since  the  area 
of  the  thin  metal  1 ization  was  increased  significantly,  increasing 
the  number  of  energetic  electrons  penetrating  into  the  depletion’ 


region.  Bond  wires  0.001  inch  in  diameter  are  used  to  contact  the 
bond  area  to  the  C3  package. 

On  the  anode  side  of  the  diode  an  AuGe-Ni  and  Au  ohmic  contact 
is  made.  This  low  resistivity  contact  is  then  soldered  to  the 
center  post  of  the  C3  diode  package.  Figure  3.8  shows  a  packaged 
diode  mounted  on  the  center  post.  Parasitic  impedances  must  be 
examined  to  understand  the  phase  characteristics  of  the  packaged 
sample. 

Figure  3.9  shows  the  lumped  element  model  for  the  measurement 
sample.  The  current  source,  IA,  is  the  induced  current  due  to 
charge  carrier  motion  through  the  sample.  The  bond  wire  inductance, 
Lb,  is  estimated  to  be  0.75  nH  for  two  parallel  bond  wires.  Package 
capacitance,  Cp,  for  the  C3  package  is  about  0.024  pF.  The  load 
impedance,  ZL,  shown  is  a  combination  of  the  fringing  capacitance  of 
the  mount  and  the  detection  circuit  input  impedance  as  seen  at  the 
plane  of  the  diode  mount.  All  of  the  passive  elements  are  constant 
with  bias  except  the  ring  capacitance.  The  inactive  region  beneath 
the  bonding  area  looks  like  a  parallel  plate  capacitor  whose 
capacitance  decreases  with  increasing  reverse  bias.  The  ring  region 
is  inactive  since  the  thick  metallization  prevents  charge  from 
penetrating  into  the  depletion  region  beneath  the  bonding  surface. 
At,  or  close,  to  punch-through  the  ring  capacitance  can  be  assumed 
to  be  constant.  Previous  measurements  attempts  were  made  at 
punch-through  so  varying  parasitic  capacitance  was  not  a  factor. 
Understanding  of  the  effect  of  packaging  parasitics,  however,  is 


Figure  3.8  Photograph  of  a  packaged  GaAs  sample 
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Figure  3.9  Lumped  element  representation  of  the 

packaged  sample  with  an  injected,  modulated 
beam. 


I 


45 

necessary  in  the  interpretation  of  the  external  current 
phase-voltage  relationship. 

Appendix  7.4  shows  the  analysis  used  to  derive  the  following 
relationship  between  the  current  source  to  be  determined  and  the 
detected  signal 


"  i-«2crVJ^VV“VrS}  (3-1) 

This  equation  shows  that  for  a  constant  ring  capacitance  (Cr)  the 
phase  of  the  impedance  multiplication  factor  is  constant.  Accurate 
measurements  must  therefore  be  done  at  or  near  punch  through  to 
minimize  circuit  related  phase  variation.  This  also  suggests  that 
calibration  with  an  empty  C3  package  would  be  futile  since  the  phase 
of  the  parasitic  circuit  would  not  be  taken  properly  into  account. 
If  the  velocity  is  assumed  to  go  to  zero  as  the  electric  field  does 
then  a  reference  point  is  achieved  for  these  measurements. 
Calibration  depending  on  other  data  is  therefore  unnecessary. 

The  packaged  diode  is  placed  in  the  sample  holder  as  shown  in 
Figure  3.10.  From  sample  capacitance  vs.  voltage  measurements  and 
from  estimations  of  the  bond  wire  inductance  it  was  determined  that 
the  source  impedance  of  the  diode  was  on  the  order  of  j5  ohms.  To 
better  match  the  diode  to  the  50  ohm  impedance  of  the  measurement 
system  a  15  ohm  quarter-wave  transformer  was  incorporated  into  the 
holder.  Teflon  dielectric  loadir?  was  used  to  reduce  the  physical 
length  of  the  transformer  and  to  act  as  a  stand  off  to  prevent 
contact  to  the  outer  conductor  wall.  The  center  conductor  of  an  SMA 
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connector  was  soldered  to  the  transformer  and  its  outer  conductor 
soldered  to  the  copper  sample  chamber.  The  detected  12  GHz  signal 
improved  greatly  with  tne  introduction  of  the  transformer. 

The  sample  chamber  mounted  to  the  vacuum  flange  is  pictured  in 
Figure  3.11.  The  diode  package  is  placed  into  the  mount  and  an 
aperture  wall  is  screwed  in  to  make  good  contact  with  the  diode 
package.  The  RF/bias  feedthrough  can  be  seen  making  contact  to  the 
back  of  the  sample  holder:  the  pipes  which  support  the  holder  are 
also  shown.  These  pipes  lead  to  a  hollow  cavity  around  the  sample 
chamber,  through  which  hot  or  cold  fluids  can  be  pumped  to  regulate 
sample  temperature.  Figure  3.12  shows  the  assembly  mounted  on  the 
far  side  of  the  vacuum  tee. 

3.4  DC  CHARATERIZATION  OF  THE  APPARATUS 


In  order  to  determine  the  performance  of  the  system  several  DC 
measurements  were  done  on  both  the  sample  and  a  blank  C3  package. 
Performance  curves  were  generated  to  help  understand  any  differences 
from  predicted  behavior  of  the  injection  or  detection  scheme.  The 
samples  used  in  these  experiments  were  GaAs  Schottky  diodes. 

Detected  current  vs.  beam  potential  is  shown  in  Figure  3.13. 
The  lower  trace  is  the  maximum  current  that  could  be  obtained 
through  the  blank  package  and  the  upper  curve  is  for  a  Schottky 
sample.  Both  have  a  peak  about  -6.5  kV  which  is  expected  since  this 
is  near  the  operating  point  and  all  the  gun  potentials  were  adjusted  - 
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Figure  3.12  Vacuum  envelope  asseml 


Sample  Current  (pA) 


Beam  Potential  (~kV) 


Figure  3.13  Detected  current  vs.  beam  potential 
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for  maximum  current  density  at  this  point.  The  ratio  of  the  sample 
current  to  the  blank  current  is  the  current  multiplication  factor. 
Figure  3.14  shows  its  experimental  relationship  to  beam  potential. 
Since  the  area  of  beam  acceptance  of  the  sample  is  less  than  half 
that  of  the  empty  package  the  peak  should  be  at  the  point  of  focus 
on  the  sample  face.  The  multiplication  factor  peak  is  therefore  due 
to  most  of  the  electron  beam  crowding  into  the  narrow  (0.02  inch 
diameter)  circle  that  is  the  sample  target.  Division  of  the  beam 
energy  by  the  band  gap  energy  would  lead  to  an  estimate  of  the 
current  multiplication  of  several  thousand. 

Reverse  bias  measurements  were  performed  on  the  sample  to 
estimate  the  limit  on  the  system  by  operating  the  Schottky  diode  at 
high  bias.  Figure  3.15  shows  the  measured  leakage  current  of  the 
sample  with  no  electron  injection.  The  upper  trace  is  the  maximum 
current  available  through  the  sample  with  a  -6.75  kV  electron  beam 
incident.  If  the  leakage  current  is  subtracted  from  the  sample 
current  the  resultant  current  due  to  injected  carrier  excitation 
looks  nearly  independent  of  the  bias.  Figure  3.16  shows  this 
result.  This  is  understandable  as  the  generation  of  electron  hole 
pairs  which  determines  the  carriers  available  is  independent  of  the 
reverse  bias.  The  increase  seen  at  -36  V  is  probably  from  avalanche 
gain  in  the  presence  of  injected  carriers. 

These  results  indicate  that  the  injection  system  and  sample 
were  operating  properly.  A  beam  current  of  330  nA  was  measured; 
however  detecting  the  RF  was  much  more  difficult. 
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Figure  3.14  Current  multiplication  vs.  beam  potential 
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3.5  TUNING  AND  AMPLIFICATION  OF  THE  DETECTED  SIGNAL 

3.5.1  Double  Stub  Tuner 

To  maximize  the  test  signal  power  detected  some  tuning  had  to 
be  provided  to  improve  the  match  between  the  source  impedance  of  the 
sample,  in  the  sample  holder,  and  the  input  to  the  first  stage 
amplifier.  A  tuner  was  needed  which  would  transmit  DC  for  sample 
bias  and  minimize  second  harmonic  pickup  from  the  6  GHz  source.  The 
detected  signal  is  about  -70  dBm  so  pickup  of  the  12  GHz  frequency 
at  this  level  will  be  noticable. 

The  tuner  used  was  a  50  ohm  airline  about  four  wavelengths  long 
at  12  GHz  with  two  25  ohm  quarter  wavelength  slugs  on  the  center 
conductor.  These  slugs  were  free  to  slide  along  the  center 
conductor  to  provide  the  proper  impedance  transformation  for  maximum 
signal  power.  Holes  were  drilled  and  tapped  in  the  slug  walls  for 
plastic  screws  which  protruded  through  a  slot  in  the  outer 
conductor.  These  screws  were  used  both  to  adjust  position  and  as 
set  screws  to  fix  the  slugs  in  place.  This  slot  in  the  tuner  body 
was  covered  with  a  piece  of  Ecosorb  with  a  slit  along  the  airline 
axis  to  accomodate  the  set  screws.  Type  N  connectors  were  used  at 
the  tuner  ends.  The  tuner  performed  adequately  with  little  unwanted 
pickup.  Machine  drawings  of  the  body,  center  conductor  and  tuning 
slugs  are  in  Appendix  2. 
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3.5.2  Cascaded  Amplifiers 


The  detected  signal  from  the  sample  diode  can  be  as  low  as  -75 
dBm.  The  HP  8411A  harmonic  converter  has  a  detection  threshold  of 
-70  dBm.  So  low  noise  amplification  was  needed  to  provide  adequate 
signal  for  phase  detection.  In  this  system  two  narrow-band,  low 
noise  amplifiers  were  used  to  raise  the  signal  to  about  -25  dBm. 

The  first  stage  is  an  11.75  to  12.25  GHz  FET  amplifier  with  28 
dB  gain  across  the  band  and  a  4  dB  noise  figure.  The  unit  was 
specially  made  with  a  bias  tee  on  the  input  for  sample  diode 
biasing.  No  separate  bias  tee  is  needed. 

The  second  stage  is  also  an  FET  amplifier  with  20  dB  gain  over 
11.25  to  12.2  GHz  with  a  noise  figure  below  6  dB.  The  two  stages 
together  provide  about  45  dB  of  gain  with  a  noise  figure  of  6  dB. 

3.6  MEASUREMENT  SCHEME 


The  measurement  circuit  shown  in  Figure  3.17  is  the  one  used  to 
determine  the  phase  vs.  applied  voltage  relationship  for  the 
injected  Schottky  diode.  The  modulation  signal  is  doubled  and  taken 
as  a  reference  signal  for  the  ANA.  The  electron  beam  is  modulated 
and  generates  charge  carriers  in  the  target  material  inducing  a  12 
GHz  current  in  the  external  circuit.  This  test  signal  is  monitored 
with  a  wavemeter  to  measure  signal  frequency  and  a  spectrum  analyzer 
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Figure  3.17  6  GHz  Measurement  circuit 
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to  observe  the  signal  level.  The  sample  is  biased  by  means  of  a 
precision  step  voltage  supply  to  increase  the  measurement  speed  and 
decrease  phase  drift  effects  due  to  beam  voltage  variation  with 
time. 

Extrapolation  to  zero  velocity  was  used  so  calibration  was  not 
necessary.  These  measurements  were  done  totally  under  computer 
control  which  generates  a  phase  vs.  voltage  table.  The  numbers  are 
then  entered  into  the  numerical  analysis  program  which  iterates 
Evans  and  Robson's  results  (3)  and  generates  a  velocity-field  curve 
which  satisfies  the  raw  data. 

Data  taken  with  this  technique  is  presented  in  Chapter  4. 
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4,  RESULTS 


4.1  MEASURED  RESULTS 


Using  the  method  and  apparatus  described  in  the  preceeding 
chapters,  the  phase  vs.  reverse  bias  relationship  was  examined  for  a 
GaAs  sample.  The  samples  used  were  Schottky  diodes,  labelled  2A26 
material,  with  a  thin  circular  metallization,  0.025  inches  in 
diameter,  and  a  small  square  bond  pad,  0.002  x  0.002  inches,  for 
contact  to  the  package  terminal.  The  capacitance  vs.  reverse  bias 
analysis  of  the  material  before  packaging  yielded  the  curve  in 
Figure  4.1.  From  this  curve  the  doping  profile  was  determined,  as 
illustrated  in  Figure  4.2.  The  epitaxial  layer  was  found  to  be  1.6 
microns  thick  with  a  doping  concentration  of  about  4x10 15  cm'3  so 
only  a  small  voltage,  about  6.5  V,  was  needed  to  punch  through  the 
epitaxially  grown  region.  Also  shown  is  the  doping  profile  for  a 
second  sample,  M13.  Figure  4.3  shows  the  calculated  internal 
electric  field  profile  for  the  reverse  bias  values  shown. 

To  test  the  Schottky  contact,  I-V  curves,  such  as  Figure  4.4, 
were  measured  to  examine  the  limits  of  reverse  bias  and  the  quality 
of  the  diodes.  Very  leaky  contacts  could  not  be  used.  The  most 
frequent  cause  of  failed  samples  was  bad  bond  wire  contacts  on  the 


Capacitance  (pF) 


Figure  4.1  Sample  2A26  capacitance  analysis. 
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Figure  4.4  Sample  2A26  current-voltage  relationship. 
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evaporated  aluminum.  Most  samples  appeared  as  open  circuits.  Other 
fabrication  difficulties  produced  rather  poor  Schottky  contacts  so 
only  two  out  of  ten  samples  made  were  usable.  The  I-V  curve  in 
Figure  4.4  was  measured  on  the  2A26  #2  sample  which  was  used  in 
these  measurements.  The  limit  of  reverse  bias  for  this  sample  was 
11  volts  as  leakage  currents  were  too  large  at  higher  biases. 

The  sample  was  placed  in  the  measurement  circuit  depicted  in 
Figure  3.17.  Phase  vs.  voltage  data  was  taken  at  300  K  with  a  beam 
potential  of  6.85  kV.  Table  4.1  shows  the  measured  phases  along 
with  their  applied  voltages.  To  establish  a  starting  point  for 
analysis,  Evans  and  Robson's  (3)  data  was  integrated  using  equation 
(2.18)  on  the  electric  field  interval  of  the  first  applied  bias. 
The  phase  which  was  generated  was  taken  as  that  of  the  first 
measurement.  All  the  other  phases  were  offset  by  a  constant  equal 
to  the  difference  between  the  measured  first  phase  and  the  generated 
first  phase.  Table  4.1  also  shows  the  adjusted  phases  which  were 
entered  into  the  numerical  analysis  program  along  with  the  applied 
voltages  and  the  sample  doping  profile.  This  adjustment  was 
necessary  to  speed  up  the  numerical  analysis  and  should  have  no 
effect  on  the  v(E)  curve  if  extrapolation  to  zero  velocity  for  zero 
electric  field  is  assumed. 

This  phase- voltage  data  along  with  the  doping  profile  were 
entered  into  the  Newton- Raphson  numerical  analysis  program  on  the  HP 
9825S  desk  top  calculator  which  augmented  the  initial  v(E)  guess  to 
fit  the  measured  data.  Figure  4.5  shows  the  curve  generated  by  the 
program.  The  flat  region  from  26  kV/cm  to  80  kV/cm  is  a  result  of  a 
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lack  of  resolution  in  the  measurements  in  this  area.  In  Figure  4.3 
the  electric  field  profile  for  the  highest  voltage,  10.4  V,  and  the 
smallest,  5.0  V,  applied  reverse  bias  are  shown.  The  lowest 
electric  field  value  for  10.4  V  is  26  kV/cm  and  the  highest  field 
intensity  for  5.0  V  is  78  kV/cm.  This  area  of  the  v(E)  curve  is 
common  to  all  the  measurements  and  changes  made  in  its  shape  will 
have  little  effect  on  the  difference  between  generated  phases.  The 
curve  in  this  area  can  only  be  offset  up  or  down  but  its  shape 
cannot  be  determined  for  normal  doping  profiles.  Reverse  bias 
measurements  were  not  made  below  5.0  V  because  of  varying  bond  area 
capacitance  and  short  depletion  lengths.  Noise  in  the  circuit 
prevented  examination  of  depletion  lengths  less  than  1.5  microns 
because  the  induced  phase  shift  is  too  small. 

The  large  discontinuities  in  the  v(E)  plot  at  certain  points 
are  probably  due  to  compensation  for  the  flat  section  of  the  curve 
in  the  numerical  analysis.  If  a  smooth  line  is  drawn  through  the 
average  of  the  plot,  a  more  reasonable  estimate  of  the  curve's  shape 
can  be  made.  The  curve  seems  to  exhibit  an  asymptotic  decrease  at 
high  fields  suggesting  that  there  is  a  saturation  velocity. 

The  amount  of  error  in  the  high  field  section  of  the  curve  may 
be  as  much  as  ten  percent.  Phase  variation  with  the  small  change  in 
velocity  in  this  region  is  difficult  to  determine.  This  phase 
uncertainty  along  with  the  resolution  of  the  ANA  (about  3  degrees) 
weakens  faith  in  this  data.  Longer  samples  would  be  useful  to 
better  the  phase  resolution  for  small  changes  in  velocity. 
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An  expansion  of  the  low  field  region  (0  to  12  kV/cm)  or  the 
v(E)  curve  is  shown  in  Figure  4.6.  The  resolution  is  better  in  this 
area  although  numerical  jumps  in  the  curve  are  still  evident.  The 
peak  velocity  is  1.5xl07  cm/s  at  3.6  kV/cm  and  the  numerical 

deviations  were  less  drastic.  The  negative  differential  mobility  is 

2 

measured  to  be  about  400  cm  /Vs  at  7.5  kV/cm. 

Finer  electric  field  increments  could  have  been  made,  however, 
limits  on  the  measurement  system's  resolution  and  on  the  number  of 
measurements  the  calculator  could  analyze  made  this  prohibitive. 


4.2  ANALYSIS  OF  WINOHORN'S  DATA 


Phase  vs.  applied  voltage  measurements  were  performed  by  Thomas 
H.  Windhorn  at  the  University  of  Illinois  (4,5).  With  his 
cooperation  his  phase- voltage  and  doping  profile  data  were  analyzed 
using  the  Newton- Raphson  iterative  technique  described  earlier. 
Figure  4.7  shows  the  electric  field  profiles  for  his  samples  at  the 
reverse  bias  values  shown.  Table  4.2  presents  the  measured  phase 
with  respected  to  applied  voltage.  The  phases  were  adjusted  as 
described  in  section  4.1  to  speed  up  the  numerical  analysis.  The 
adjusted  phases  in  Table  4.2  along  with  the  applied  voltages  and  the 
doping  profile  were  entered  into  the  numerical  analysis  program  and 
the  result  was  the  curve  in  Figure  4.8.  Again,  numerical  deviations 
are  evident.  However,  the  field  variation  across  the  sample  is 
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Figure  4.6  Low  field  velocity  curve 
Close  up  of  Figure  4.5. 
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Table  4.2  Adjusted  phase  vs.  applied  voltage 
from  T.  Windhorn's  measurements. 


V 

4>  adjusted 

<t>  measured 

-1.0 

-  .713 

.593 

-2.0 

-  .826 

.480 

-3.0 

-  .939 

.367 

-4.0 

-1.044 

.262 

-5.0 

-1.114 

.192 

-6.0 

-1.184 

.122 

-7.0 

-1.227 

.079 

-8.0 

-1.262 

.044 

-9.0 

-1.219 

.017 

-10.0 

-1.306 

-.00001 

-15.0 

-1.393 

-.087 

-20.0 

-1.463 

-.157 

-25.0 

-1.511 

-.205 

-30.0 

-1.559 

-.253 

-35.0 

-1.598 

-.202 

-40.0 

-1.646 

-.340 

-45.0 


-1.699 


-.393 


0*7  cm/s) 
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small  enough  (20  kV/cm)  and  the  Schottky  was  very  good  so  high 
reverse  biasing  could  be  achieved.  Therefore,  electric  field  limits 
were  calculated  across  the  entire  area  of  numerical  analysis.  No 
large  flat  region  is  observed  and  an  average  line  is  in  closer 
proximity  with  the  actual  curve.  The  smooth  curve  agrees  with  Dr. 
Windhorn's  published  results  exactly. 


4.3  COMPARISON  OF  OUR  RESULTS  WITH  WINDHORN'S  RESULTS 


The  purpose  of  T.  Windhorn's  data  was  to  measure  the  high  field 
region  of  the  v(E)  curve.  These  measurements  are  too  coarse  to  be 
used  for  examination  of  the  curve  at  or  about  the  peak  velocity 
point.  Comparisons  should,  therefore,  only  be  made  with  other  data 
at  high  fields  (>10  kV/cm). 

Figure  4.9  is  an  overlay  of  the  v(E)  curve  as  generated  from 
Windhorn's  phase-voltage  data  and  that  which  satisfies  our  data.  At 
electric  field  values  within  20  kV/cm  of  the  peak  the  curves  track 
within  less  than  10  percent  of  each  other.  At  higher  fields  (>50 
kV/cm)  Windhorn's  data  indicates  a  constantly  decreasing  velocity 

with  electric  field.  Our  curve  begins  to  flatten  and  approaches  a 

6 

constant  velocity  of  about  8.8x10  cm/s  above  80  kV/cm.  This  is  a 
major  discrepency  as  to  the  value  of  the  second  valley  mobility  and 
where  it  may  become  noticeable.  Further  experimentation  is 


Velocity  (10*7  cm/m) 
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required  to  understand  this  difference  in  the  results  or  to  prove 
which  relationship  is  correct. 
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5. CONCLUSIONS 


5.1  LIMITATIONS  OF  THE  UNIFORM  FIELD  APPROXIMATION 


Although  it  has  been  shown  that  this  is  not  the  case  for  the 
samples  used  in  these  experiments,  a  velocity-field  relationship  can 
be  derived  from  the  ..base- voltage  data  in  Table  4.1  using  the 
uniform  field  approximation.  The  electric  field  is  taken  as  the 
applied  bias  voltage  divided  by  the  epitaxial  layer  length  and  the 
velocity  is  calculated  using  equation  1.6.  Figure  5.1  shows  the 
resulting  uniform  field  curve  and  the  curve  generated  using  a 
non-uniform  field  numerical  method. 

The  electric  field  variation  across  the  depletion  region  in  the 
punched- through  case  for  our  sample  was  about  90  kV/cm.  The  uniform 
field  approximation  tends  to  compress  the  range  of  measurement  in 
the  electric  field  and  arrive  at  an  average  velocity  for  the 
interval  of  the  v(E)  curve  being  sampled  in  the  measurement  point. 
In  comparison  the  uniform  curve  demonstrates  a  higher  velocity 
compared  to  the  same  interval  on  the  non-uniform  curve.  This 
demonstrates  the  increased  accuracy  of  the  non-uniform  field 
numerical  method  over  the  simplified  uniform  field  approximation. 
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E lac trie  Field  (kV/c*) 


Figure  5.1  Comparison  of  non-uniform  field  solution 
with  uniform  field  solution. 


5.2  SUMMARY 


To  measure  the  velocity  of  electrons  drifting  through  a 
semiconductor  material  relative  to  the  internal  electric  field,  a 
microwave  frequency  modulated  electron  beam  injection  system  was 

built  and  demonstrated.  In  order  to  permit  the  examination  of 
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samples  with  high  doping  concentrations  (>10  cm  ),  the  effect  of 
a  non-uniform  electric  field  taper  was  accounted  for  in  the  analysis 
of  the  measured  signal  phase-voltage  data.  The  numerical  techniques 
employed  required  shorter  samples  than  those  usually  used  in 
previous  v(E)  measurements.  For  adequate  phase  resolution  high 
carrier  injection  rates  were  necessary  (12  GHz). 

An  electron  gun  and  reliable  high  voltage  supply  were  assembled 
to  provide  a  330  nA  beam.  An  efficient  6  GHz  modulator  was  designed 
and  employed,  giving  good  single  frequency,  single  mode  transverse 
beam  modulation.  RF  power  was  transmitted  through  the  glass  vacuum 
envelope  so  the  vacuum  was  not  compromised.  Amplification  of  the 
low  power  12  GHz  test  signal  was  necessary  to  measure  the  injected 
sample  current  phase-reverse  applied  voltage  relationship. 

An  examination  of  the  sample  package  paras i tics  was  performed 
to  determine  non-punched-through  phase  behavior.  The  voltage 
dependance  of  the  bond  pad  area  capacitance  (Cr)  and  the  size  of  the 
bond-wire  inductance  were  discovered  to  be  non-negligable. 
Calibration  with  an  open  package  would  not  take  these  impedances 
into  account  and  would  therefore  be  useless.  A  zero  velocity  was 
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assumed  for  a  zero  applied  electric  field  so  a  reference  was 
established  for  the  v(E)  curve. 

Implementation  of  this  scheme  yielded  a  good  estimation  of  the 
v(E)  curve  around  the  peak.  This  data  correlates  well  with  other 
low  field  measurements.  Numerical  analysis  of  T.  Windhorn's  (4,5) 
data  has  demonstrated  good  agreement  with  the  uniform  field 
approximation  for  lightly  doped  material  at  high  electric  fields. 
The  electric  field  variation  across  the  sample  length  was  small  ( <20 
kV/cm)  so  for  areas  on  the  v(E)  curve  of  small  velocity  variation 
with  electric  field  the  approximation  is  appropriate. 

Good  samples  were  not  available  in  quantity  before  the  end  of 
this  project.  Truely  conclusive  data,  therefore,  could  not  be 
obtained  about  the  high  field  nature  of  the  generated  v(E)  curve. 
Further  experimentation  is  needed  to  establish  the  behavior  of  the 
drift  mobility  for  high  fields. 

Examination  of  the  electric  field  profiles  of  higher  doped 
semiconductor  samples  has  led  to  the  design  of  samples  to  improve 
the  range  and  accuracy  of  measuring  the  v(E)  relationship  using  the 
numerical  analysis  method  presented  here.  To  get  high  enough 
resolution  for  close  examination  of  the  peak  of  the  v(E)  curve, 
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lightly  doped  (about  2x10  cm  )  material  is  necessary  to  decrease 
the  electric  field  taper.  Epitaxial  layer  lengths  of  about  5  ^m 
would  provide  a  field  taper  of  about  10  kV/cm  to  permit  detailed 
estimation  of  the  low  field  influence  on  electron  velocity.  Smaller 
bond  pad  areas  will  help  decrease  the  influence  of  parasitic 
capacitance  variation  on  the  phase.  Experience  in  making  these 
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measurements  has  led  to  a  few  specific  suggestions  to  improve  the 
technique  for  higher  doped  material. 


5.3  RECOMMENDATIONS 


A  more  lightly  doped  sample  (about  7xl014  cm"3)  with  a  longer 
epitaxial  layer  (about  4  microns)  would  allow  enough  field  variation 
to  examine  the  v(E)  curve  up  to  about  200  kV/cm.  The  longer  sample 
lengths  would  also  allow  larger  phase  shifts  for  smaller  electric 
field  increments  increasing  the  resolution  of  the  resulting  curve. 
Very  good  Schottky  contacts  are  also  necessary  for  high  reverse 
biases  to  improve  the  range  of  measurement. 

Numerical  analysis  should  be  executed  on  a  more  powerful 
computer  than  the  HP  9825.  Ideally  an  intellegent  choice  of  bias 
points  can  minimize  the  number  of  measurements  by  increasing  the 
resolution  around  regions  of  high  mobility  variation  with  electric 
field  (i.e.  the  knee  of  the  curve  for  III-V  materials)  and  having 
lower  resolution  in  areas  of  less  velocity  variation.  Even  with 
this  selective  method  of  examination,  the  number  of  measurements 
required  to  analyze  a  large  part  of  the  curve  would  be  about  200. 
With  specialized  numerical  analysis  techniques  and  advanced  program 
chaining  methods  the  memory  requirements  and  sheer  computing  power 
necessitates  the  use  of  a  minicomputer  if  not  a  main  frame. 
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Diffusion  characteristics  for  semiconductors  is  also  important 
in  the  understanding  of  device  behavior.  To  measure  its  dependence 
a  means  of  varying  the  charge  pulse  width  and  separation  within  the 
depletion  region  would  be  ideal.  A  modulator  with  a  large  frequency 
range,  say  6  to  8  GHz,  would  be  useful  to  vary  the  pulse  injection 
rate.  The  variation  in  charge  separation  along  with  knowledge  of 
the  v(E)  curve  should  yield  enough  information  to  estimate  a 
diffusion-electric  field  relationship. 

With  these  and  other  appropriate  improvements  the  CW 
time-of-flight  measurement  technique  should  be  applicable  to  a 
variety  of  measurements  on  a  any  semiconductor  material  on  which  a 
non-conducting  contact  can  be  made  and  an  internal  electric  field 
can  be  established. 
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0.  * 11863C  SOFTWARE  FOR  SYSTEM  FOR  TRANSPORT*! 
li  'ANALYSIS  OF  SEMICONDUCTORS,  4/8/82* i 

2i  *P/N  11863-10003,  MASH I NGTON  UNIVERSITY  MICROMAVE  LAB* • 

3>  ‘MODIFIED  BY  OLEH  HOMSEY  ON  DISK  UNDER  NeaTaa*. 

4i 

Si  dap  *TAS  MEASUREMENT  PROCRAM* (Milt  2000 
6a 

7 •  ‘TEST  CONDITIONS  SECT  I » , 
8.  'INITIALIZATION* ■ 

9i 

10.  ant  *Entar  tha  nuaber  of  aaaa.  pta.'.H 
lli  dla  X(H),YtN],EI2),U(N],ZtH) 

12.  dla  AI61, 8161, D(N),KI2I. 0131, OINI ,TI31 ,VIN) 

13.  dla  E«(7],P*I4],S«I20] 

14. 

15.  dav  * A* ,710 , *S* ,706, *R* ,716, *M* ,714, *6* ,704 , *P* , 713 

16. 

17.  ‘CAPACITANCE  OF  OPEN* ■ . 0817KI 1 1 

18. 

19.  ‘SYSTEM  IMPEDANCE*. 301KC2I 

20. 

21.  radiell  7|lcl  7jr«a  7 

22. 

23.  ‘CALIBRATION  SECT  I  ONnannanaaan—anneanannananaannnne—nneenennannnneaan* 

24. 

23 >  beepicfg  3, dap  ‘SET  TEST  CHANNEL  GAIN  TO  60  dB'iatp 

26.  ‘STORE  BEAM  CENTER  OFFSETS*. 

27.  baaptdap  ‘HOLD  BEAM  CTR  and  PRESS  CONTINUE* iatp 

28.  DOlZigab  ‘Avar* 

29.  V10I11 i2>Zigab  ’Avar* 

30.  V)0(2] 

31.  baapidap  *Plaea  10V  Rev.  8laa  on  Dlode'iatp 

32.  DBigab  ‘Attan* 

33.  cfg  7)1)1 igab  'Read* 

34.  Oat n*C. 5(6-4) )/B)0 131 iafg  7 

33.  0)0 

36.  DEC  11  |0>EC21  igto  *Maaa* 

37.  ant  'la  Callb.  Data  on  Tape?  CY/N)*,P* 

38.  If  eapCP«Il))-*Y*igto  'Caload* 

39. 

40.  *S21eal* • 

41.  gab  'Thru* 

42.  ant  *Mant  To  Store  Callb.  Data?  <Y/N7*,P» 

43.  If  cap(P*Cl))«*Y*igto  'Calator* 

44. 

43.  ‘MEASUREMENT  SECT  ; OMananaaaaaaaaaaaaaaaaaaaaaaaaaaaaamaannaaannaaaaaaa* , 

46. 

47.  *Maae*. 

48.  DC|l>6(baapiant  'Enter  Saaple  Length  Cua)*,H 

49.  Hale-6>H 

50.  baapidap  'Set  Up  Device  Under  Teat'iatp 

31.  ant  ‘Enter  Modulation  Fraq  CMHi)‘,F 

32.  ‘MaaaAgn*. 

S3,  afg  3 

34.  beepiant  ‘Enter  Blaa  CVol tal * ,V(CI 
S3,  F1ZIC1 

36.  8>I 

37.  gab  ‘Loop* 

38.  ell  *Rtop'<XtC],YCCl,D(Cl,Q{Cl) 

39.  If  C>Ntgto  ‘Table* 


C*l>Cigto  ’MeasAgn 


60 
61 

62 1  ‘MAIN  MEASUREMENT  LOOP— — H— 

63. 

64.  ‘Loop*. 

65.  gsb  'Reed* 

66.  If  I>6i99b  ’Noraal* 

67 .  rat 

68. 

69.  ’CHECK  DATA  SECT  1  OWnamni 

70. 

71.  ’Check* > 

72. 

73.  ell  'Rtop'IXlCJ,YlC1,0IC1,atCl> 

74.  -<B2c6MHbZ(C]/0[C]}U(C1|flt  3 

75.  dsp  ’VELOCITY  •• ,UCC1 , *e/sec * ibeepi stp 

76.  ant  ’Meat  to  retake  this  pt?  (Y/H>*,P* 

77.  If  eaptPOf 1 J)**Y* 19 to  ’MeaaAgn* 

78.  If  C>N,gto  ’Table* 

79.  ant  ’Continue  to  next  pt?  CY/NJ’,P* 

80.  If  cep<P0lll}«*Y* |C»1 >Cigto  ’MeasAgn’ 

81.  dep  ’MAITIHG*  »xtp 

82.  C’llC.gto  ’MeasAgn* 

83. 

84.  ’DATA  OUTPUT  SECTIOMaaBBaain— Baanaasanaseni 

85. 

88.  ‘Table’. 

87.  ent  ’Enter  Table  Title  (<20  Char.)’, SO 

88.  fat  l,c,2/iwt  ’M.l’.SO 

89.  fat  1,’OATA  PT’.ziart  ’M.l* 

90.  fat  2, 8x. ’PHASE  Crad)* ,Z|*rt  ’M.2* 

91.  fat  3,4x, 'MAGNITUDE', ziart  'W.3* 

92.  fat  4,6x,’FREQ  (MHz) * , Z i«rt  *M.4> 

93.  fat  5,6z, ’BIAS  <V>’,Zi«rt  *M.S* 

94.  fat  6, 7x, ’LENGTH  (ual’iwrt  ’M.6* 

95.  for  C“1  to  N  by  1 

96.  fat  l,3x,f2.0,zi«rt  ’M.l’.C 

97.  fat  2,10z,f6.4,zi«rt  *M.2’,0tCl 

98.  fat  3,10x.fS.3,zi*rt  *M.3*,0IC] 

99.  fat  4,10x,f5.0,zi«rt  *M.4’,ZtCI 

100.  fat  5,10z,f7.3,zi»rt  *M.5*,VtC] 

101.  fat  6,9x,f9.2|«rt  *M.6’,Hale8 

102.  next  C 

103. 

104.  dep  ’DOME*  i»tp 

105. 

106.  ’SUBROUTINES  SECT  I  ONsaeassssasas—nsnaaw 

107.  ’Caload* ■ 

108.  beepidnp  ’Insert  Data  Tape’iwalt  2000 

109.  ent  ’Enter  Track  to  Find  Cel lb.  Data*, I 

110.  ent  ’Enter  File  to  Find  Callb.  Data*,J 

111.  trk  1 1 ldf  J.EIel 

112.  9to  ’Meae* 

113. 

114.  ’Calstor*. 

115.  beepidnp  ’Insert  Data  Tapa’ivait  2000 

116.  ent  ’Enter  Track  To  Store  Cel  lb.  Data*, I 

117.  ent  ’Enter  File  To  Store  Callb.  Deta’.J 

118.  trk  tiref  J.ECSI 

119.  gto  ’Maas’ 

120. 

121.  * 

122.  ’ 


Reed*. 

Chi’  1 llZ.gsb  ’Aver 
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123.  V-OUDE 

124.  *Ch2* .2)Zig9b  ’Avar* 

125.  ell  *Rtop' CE.V-0E2I >9, A) 

126.  If  8>2.2igto  *Up* 

127.  If  8<.4igto  "Down" 

128.  *Ch3* < If  0<3  and  f l97|0I3]B/tn»<. 5(6-6)) >Digto  *3 

129.  4>Zi9sb  ‘Avar* 

130.  tn»CV-. 5(6-0)) >0 

131.  ell  'Ptor'CO.A.X.Y) 

132.  If  fIgSigto  ‘Data* 

133.  *Cal*.X>EEll(V>E(2)irat 

134.  ‘Data* .X1XEC1 |Y>YEC1 .rat 

135.  ‘Up*. 

136.  0*1>G| If  Q>6i6>0iprt  * OVERRANGE* »gto  *Ch3* 

137.  99b  *Attan* 

138.  gto  'Raad* 

139.  ‘Down*. 

140.  8-l>6«lf  0<l|l>0i9to  *Ch3* 

141.  99b  ‘At ton* 

142.  gto  ‘Rood* 

143. 

144.  'AVERAGE  MEASUREMENTS* ■ 

145. 

146.  ‘Avar*. 

147.  OMItSOOH) 

148.  for  L*1  to  D 

149.  99b  ‘Volt* 

150.  U*V>U 

151.  naxt  L 

152.  (U/D).0025>V 

153.  rot 

154. 

155.  ‘59313A  A/0  REAOIMO* ■ 

156. 

157.  'Volt*. 

158.  fat  l,*H*.f.0,*AJ*iart  *A.1*,Z 

159.  rot (rdb(*A*) ,8)*rdb(*A*) )V 

160.  rot 

161. 

162.  ‘SET  8418  IF  ATTENUATION*. 

163. 

164.  * At  ton* • 

165.  fxd  0 

166.  boap 1  dip  'SET  TEST  CHANNEL  GAIN  TO* , 10C7-6J , stp 

167.  rot 

168. 

169.  ‘CALIBRATION  WITH  THRU-SINGLE  PATH*. 

170. 

171.  ‘Thru*. 

172.  if  041|baopid9p  ‘Conn  Caltb  Soctlon'.stp 

173.  1H  1  l>0|9fg  lO.ggb  ‘Loop* 

174.  efg  10 

175.  rot 

176. 

177.  ‘MATH  DURING  MEASUREMENT*. 

178. 

179.  'Noraal *  1 

180.  ell  'Cdiv'(X(Cl.YtC],E(l),EE21,TE2],T[3]> 

181.  TE2]>X(C]tTI31>YEC1 

182.  rat 

183. 

184.  ‘POLAR  TO  RECTANGULAR*. 

185. 
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186.  *Ptor  *  < 

187.  plcoa(p2)>p3 
188 1  plsin(p2)>p4 
189t  ret 

190  • 

191.  * RECTANGULAR  TO  POLAR*. 

192. 

193.  *Rtop* i 

194.  \ Cplpl*p2p2) >p3isfg  14 

195.  2etn(p2/ Cpl«p3)) >p4 

196.  cfg  14 i ret 

197. 

198.  ‘COMPLEX  DIVIDE* • 

199. 

200.  *Cdlv*. 

201.  p3p3+p4p4  >p7 

202.  Cplp3«p2p4)/p7>p5 

203.  Cp2p3-plp4)/p7>p6 

204.  ret 

205. 

206.  end 
■12600 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  8uNE*U  OF  STANOAftOS  -  *963  -  A 


0 •  *v<E)  Nuaerlcal  Analysis  Prograa  for  Arbitrary  #(V)‘i 
li  ‘Developed  by  Olan  Hoasey,  7/8/82.  On  Disk  Undar  CorRaf ‘ i 
2 1 

3.  rad 
4| 

Si  dap  *v(E)  NuaAnl  Prog  for  Any  8(V)'i«eit  2000 

6i  ant  *Entar  Modulation  Frag.  In  Hi*,M 

7 1  fit  Siprt  ‘Mod  Fraq>*.M«MM2aOM 

8.  ant  ‘Entar  Huaber  of  Haas.  Pts.'.M 

9i  ant  ‘Entar  Nuabar  of  Doping  Pts.‘.H 

10 ■  ant  ‘Entar  Epllayar  Depth  In  alerona'.L 

Hi  L>le-6>l. 

12i  ant  ‘Enter  No.  of  Initial  Duass  Points*, 0 

13i  dla  C(0] ,D(H] ,E(21,N1,H(21],L(N],P(H] ,0(N1,SC01 .VtOiN).X(M) 
14 i  dla  ltNl.JCN.Nl ,VtN] ,  (1(01, F(01,G(0], 0(01,  M(01 
15.  dla  AI2.H1, 814.01 

16i  1 . 6021892s- 1910 i ant  ‘Enter  Ral.  Dial.  Const. *,r2 
17i  fxd  2iprt  ‘Ral  Dlal-‘,r2,‘  ‘ 

18.  8/(r2/(3SaOala-9)>8 
19 1 

20i  beep.dsp  ‘Doping  Data  Entry* .volt  2000 

21 i  dsp  ‘Enter  X  in  alerons  and  0  In  ce-3*inalt  2000 

22i  prt  ‘Doping  Data*,'  ■ 

23.  for  1>1  to  M 

24.  fxd  2 

23.  ant  XUliprt  l.‘X-‘.XtIl 

28.  fit  3 

27.  ant  DCUiprt  ‘0-‘,0(l],‘  • 

28.  XII ]ala~6>X( 11(0(1 ]a8aia6>0( I ] 

29.  next  1 

30.  prt  •  • 

31.  fxd  2.01F(1] 

32.  for  1*1  to  N-2 

33.  ell  ‘CorrTrap'Cl, XIII, XI 1*11 ,Xf 1*21 . 0(11, 011*11, 0(I»2J,AJ 

34.  A»F( 1 ])F( 1*11 
35 •  next  I 

36.  .Sa<D(M]*0(M-l])a(X(MJ-X(M-l])*F(M-l]>FCM] 

37 .  for  I “1  to  H-l 

38.  ell  'Peel,  in*  (X( I1.XCI»1J, Fill, Ftl-»1J,  AC  1,11,  At  2, 111 

39.  next  I 

40.  0>V 

41 i  for  K*1  to  21 

42.  YM./201Y 

43.  for  1*1  to  M-l 

44.  If  X(U<Y  and  Y<-XI  I«ll*le-9,D£  1 1»<-13>HIKJ  ,M-1H 

45.  next  I 

46.  next  K 

47. 

49.  baap.dsp  ‘Phase  vs.  Voltage  Entry*.valt  2000 

49.  dsp  ‘Entar  Phase  In  radians* .wait  2000 

50.  prt  ‘Phase* ,* Inforaat Ion* , *  * 

51.  for  I«1  to  M 

52.  ant  8(1] .prt  !,*«-•, 0(11 

53.  ant  Vdliprt  ‘V-‘,V(I1,‘  • 

34.  vm-.7smu 

35.  next  t 
38.  prt  ■  • 

37.  l>l)Jt0>X>V(01>A>Uiefg  1 

38.  ‘Step*. 

59.  X*L/300>X 
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60 1  if  X>L|9to  • Punch* 

61 .  *Stap* i 

62 ■  If  X[J«ll<X|J*l>JiV«V(I-l]>U|sfg  l.gto  'Step* 

63 i  -Xa<A[ 1 , J] *A(2 , JI>X) >F 
6ft  tf  not  flgltjap  2 

6St  (ACl,Jl*At2,J]>(X*XUl)/23B(X-XUl)>D 
66 i  If  flgl,jap  2 

67.  (A(l,J]»At2,J]M(A*X)/2)«(X-A)>0 

68.  -V(I-11«0»U>V 

69.  if  F*V<Vt I 1 icfg  1» I»l> I .F/X1EC1. 1-11 |X>A|0>Uig9b  *EGen* 

70.  If  I>N*ligto  ‘Initial* 

71.  gto  'Step* 

72. 

73.  ‘Punch*. 

74.  (Vt I 1-V)/L>£11 , I ] I I«l> I igsb  *EGen* 

73.  if  I4N»ligto  ’Punch* 

76. 

77.  *(nitial*. 

78.  baep.ont  'Enter  Data  Flle*,J 

79.  trk  O.ldf  J.Ulal.Fla] 

80. 

81.  for  1*1  to  H 

82.  Cll  'Generate' C I } 

83.  Ltmitllifzd  f.prt  Ptll 
84 ■  next  I 

83.  prt  *  * 

86.  ’Again*. 

87.  for  K-l  to  N 

88.  cll  'ttarata'(K) 

89.  for  J“1  to  N 

90.  efg  7 

91.  if  ECZt.JKZ  or  7»next  J 

92.  if  flg7|gto  96 

93.  cll  'Generate' CJJ 

94.  (LT JI- I [ J])/le3>JI J,K) 

95.  next  J 

96.  cll  'Dei tore to '<KJ 

97.  next  X 

98.  Inv  J>J 

99.  nepy  -lain 

100.  eat  AUY 

101.  for  I»1  to  N 

102.  for  S»1  to  0 

103.  if  ECl,ll<«FCSJ  and  F[S]<-Et21 , IJ »Yt IXUtSlHJIS] 

104.  next  S 
103.  next  I 

106.  1»X 

107.  for  I«1  to  N 

108.  cll  'Generate'CI) 

109.  0>Y»  if  abn(L(  I IX.  001 1 1>Y 

110.  fxd  4, prt  P(II 

111.  X«Y)X 

112.  LUlmil 

113.  next  I 

114.  prt  •  • 

115.  if  X .gto  *v<E>  Plot* 

116.  gto  ‘Again* 

117. 

118.  ’SUBROUTINES*. 

119. 

120.  'Generate* if  It  3icfg  8 

121.  Ell >pll>G(11 i 2 IT 

122.  for  S-0  to  1  by  -1 


123.  If  Etl,pll<*FCSl  and  FCS1<*E[21 ,pll |UCS1*C-IJ >MtTl 'FtSHGCTI ,afg  8 

124.  if  Etl ,plI<*FtS]  and  FtS1<*£t21,pll|T*l»Tigab  ’Slopa* 

128.  If  flg8  and  FCSl>Et21,plJ,FtSnGlT«ll,tKSl“C-l>>*«T*ll 
128.  If  l*S|]ap  2 

127.  If  f  lg8  and  FtS]>Et21,pll  »FIS-ll>GtT*2J  |UtS-ll*C-l>H«T*21  »efs  8»1>S 

128.  naxt  S 

129.  H(21>H(l]«0C21>0tl] 

130.  If  f lg8|T-l>R>T|GtT]>G(T*l] |MtT]}M(T«ll i jap  2 

131.  E(21.pl]>GtT)|CMtT-l]*MCT«lJ)/2>MCT)i0tT-ll>0(TI>0(T*l]|T»l>R 

132.  0>Eil>Ctl]|0>Stl] 

133.  far  S-l  to  R-l 

134.  ell  'CarrTrap'<pl,6(SI,GfS«ll,GtS*21,MtSl,MfS*ll,H(S*2],A) 

135.  M01S1*E>E 

136.  eog(MO>CtS»l]|9in(NME)>StS*l] 

137.  naxt  S 

138.  0>E»F 

139.  for  S-l  to  R-l 

140.  ell  'CorrTrap'<pl,GtSl,GtS*ll,GtS»21,CCS],C[S*l],CtS»2],A> 

141.  MOISUEIE 

142.  ell  'CorrTrap'Cpl,GtSl.G[S*l],GtS»2].S[S],S(S*l],StS*2],B> 

143.  8WUS1*F>F 

144.  next  S 

143 •  \ CF*2*E*2> > A i 2a t n<-F/ (E«A> ) >Ptpl ] 

148.  (aipii-PCpin/aipiiiLtpii 

147.  rot 

148. 

149.  ■PeeLln'i 

150.  (p4-p33 / (p2-pl) >p6 

151.  p3-pl9p6>p5 

152.  rot 

153. 

134 .  'CorrTrap* . 

135.  p3-p2)rliif  rl-0ir«t 

156.  Cp6-p33/rl>r2|if  p4-p3*0irl/2a(p6«p5)>p8irat 

157.  Cp7-p6)/<p4-p3)>r3 

158.  rt/2aCp6*p5)»rl*2/12«Cr2-r33>p8 

159.  rat 

160. 

161.  Mntagrata*. 

162.  p6ap2*p7xp2*2»p8ap2“3«p994>2*4}rl 

163.  p6»p3*p7ap3“2*p8ap3*3*p96p3“41r2 

164.  r2-rl»p4>p5 

165.  rat 

166. 

167.  *EGan* . 

168.  L/20>Y|l>rl|I-l>r2|2>r3 

169.  ‘Eatep*. 

170.  If  Xtrl*l]<Yirl«l>rli9to  •Estep* 

171.  Etl,r2]*A.Cl,rl]»A(2,rl]»Y>Etr3,r2] 

172.  tf  E(r3,r21>0|0>Etr3,r2Iiret 

173.  if  r3*21.rat 

174.  r3*l>r3|Y»L/20>Yigto  *Eatap« 

175. 

176.  *Itarata*  > 

177.  afg  4 

178.  far  1*1  ta  0 

179.  efg  5 

180.  if  £(l,pll<*F(n  and  Ft  Il<*Et21  ,pll  »Ul  I  J*la3»UU  1  .FtlDZ.afg  5 

181.  if  f lg4  and  f IgS.FC I llY.efg  4 

182.  if  not  flg4  and  not  flg5iO>I 

183.  if  Etl  ,pll>FC  1 1 10>1 

184.  naxt  I 

185.  rat 


186 1 

187i  'Del  terete*. 

188 i  ifj  4 

189 1  for  1*1  to  0 

190 •  of 9  5 

191 1  If  Etl.plK-FCIl  and  Ft  I  ]<*£(21,pll  .lit  1 1  -  la3>Ul  I  ]  .efg  4,cfg  S 
192 ■  if  not  flg4  and  flg5.0>I 
193 i  If  Etl,pl]>FtlliO>! 

194 i  nazt  I 
195 •  pat 
196 1 

197 >  *SIopa* • 

198.  for  U*20  to  1  by  -1 

199.  If  EtU,plJ<-FtSJ  and  FCS]<*E(U«l,pl] ,HtUl>0tT-U ,1>0 

200.  nazt  U 

201.  rat 

202. 

203.  *v<E)  Print'. 

204.  fit  3 

205.  for  1*1  to  0 

206.  If  Etl.NK-FUl  and  Et21  ,ll>*Ft  1)  iprt  'E'.FIIJ,  *V'  ,Ut  II,'  * 

207.  nazt  1 

208.  gto  *w<E>  Doap* 

209. 

210.  *vCE>  Plot'. 

211.  art  713, 'pgl.ln.lplOOO, 2000, 7200, 10000.' 

212.  -lntCEtl,Nl>la-6-l)al0)2|l>V 

213.  If  Z>20  »5>Y 

214.  if  Z750.101Y 

213.  art  713, 'ao0,' ,Z, • ,0,3. ' 

216.  art  713, *zo0 , ■ ,Y, * ,0 , • ,2, • , • 

217.  fzd  0 

218.  for  K-0  to  Z  by  Y 

219.  art  713,'po*,atr<K) ,',0.' 

220.  atb  713,*cpO,-l|lo6|lb*,ztrCK} ,3 

221 .  nazt  K 

222.  art  713, *po> ,Z/2,* ,0icp0,-3i' 

223.  atb  713, 'IbElaetrle  Field  In  KV/ca',3 

224.  fzd  2 

225.  3>A. If  0>1.371a9|S>A 

228.  art  713,'ye0..23,0,',6,'.' 

227.  for  K-0  to  •  by  .25 

228.  art  713,'paO.' .atrCKJ , ' .* 

229.  atb  713, *ep-3.5,0| loS.lb* ,atr <K) ,3 

230.  nazt  K 

231.  fzd  6 

232.  art  713, 'paO, ■ ,0/2, • , • 

233.  art  713,*cp-10,0|dl0.1i* 

234.  atb  713,'lbVeloolty  z  10*7  ea/a',3 

235.  fzd  6 

236.  for  1*1  to  0 

237.  If  Etl,Nl<-FtI]  and  EI21 ,11)*FU]  .gab  'Plot* 

238.  nazt  I 

239.  gto  *v<E)  Doap* 

240.  'Plot*. 

241.  If  not  flglO.art  713, 'papa*, -Ftllble-S, ',' ,UtIlala-5, • , • »afg  10 

242.  if  flglO.art  713,'popd'  ,-Ft  Ilble-5,  • ,  •  ,Utn«le-5,  ■ .  • 

243.  fit  3.prt  UII) 

244.  ret 

245. 

246.  *v<E)  Ouap*. 

247.  ant  'Enter  Track  For  w(E)  Storage*, J 
249.  ant  'Enter  File  For  v<E>  Storage*, K 
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249 «  trk  J,rcf  K,UO),Ftol 
2S0 i  tad 
251 1 

252 ,  'EOup'i 
2S3i  for  1-1  to  N 

254,  fxd  2iprt  V(I) 

255,  for  J* 1  to  21 

256,  fit  3,prt  EtJ.ll 

257,  if  El J, I }»0,21>J 

258,  noxt  J 

259,  prt  •  • 

260,  noxt  I 

261,  dop  *Duapcd*utp 
■17309 
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Figure  A2.2  6  GHz  Waveguide  Launcher 
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APPENDIX  3 

Power  Supply  Circuits 
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figure A2. 7  High  voltage  circuit 
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Figure  A2.9  Regulating  circuit  H.V.Il 
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APPENDIX  4 

Lumped  Element  Sample  Model 

At  the  node  called  the  package  post  point  in  Figure  3.9,  one 
can  relate  the  three  currents  through  equation  A. 4.1:  ^  is  the 
current  through  the  ring  capacitance  into  the  node;  I 2  is  the 
current  through  the  parallel  combination  of  the  package  capacitance 
and  the  load  impedance;  and  I A  is  the  excitation  current  to  be 
examined.  So 


Xl  -  *2  -  XA 


(A4.1) 


1  +  jwC  ZT 
pL 


(A4.2) 


By  applying  Kirchoff's  voltage  law  to  the  external  loop  one  finds 


3^7+Z*  *  l2  +  l»  hl2m0 


(A4.3) 


substituting  (A4.1)  for  ^yields 


(A4.4) 


1  +  J^Z'  - 


1 


The  measured  signal  is  then  just 


V»  “  Z’  *  *2 


ZLXA 


,  2  - 5 - :  (A4.5) 

1"“  CpCrV 


which  can  be  described  as 


Vm  -  *  !A  (M'6) 

Equations  A4.5  and  A4.6  demonstrate  the  significance  of  the 
bond  area  capacitance  phase  contributions  to  the  measured  phase  of 
the  signal.  As  the  applied  reverse  bias  decreases  below  the 
punched- through  value,  Cr  increases  and  changes  the  phase 
contribution  of  the  effective  impedance  in  A4.6.  To  minimize  the 
effect  of  this  capacitance  shift  the  bond  area  can  be  reduced, 
increasing  the  shunt  impedance  across  the  actual  current  source. 
Measurements  could  also  be  performed  with  reverse  applied  bias 
voltages  near  or  above  the  punched-through  value.  The  variation  of 
bond  area  capacitance  should  be  very  small  in  this  case  as  the 
depletion  depth  does  not  change  a  great  deal. 
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